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The physics impact of a staged approach for double-/? decay experiments based on ,f ’Ge is studied. 

The scenario considered relies on realistic time schedules envisioned by the Gerda and the Majo- 
rana collaborations, which are jointly working towards the realization of a future larger scale ,6 Ge 
experiment. Intermediate stages of the experiments are conceived to perform quasi background-free 
measurements, and different data sets can be reliably combined to maximize the physics outcome. 

The sensitivity for such a global analysis is presented, with focus on how neutrino flavor models can 
be probed already with preliminary phases of the experiments. The synergy between theory and 
experiment yields strong benefits for both sides: the model predictions can be used to sensibly plan 
the experimental stages, and results from intermediate stages can be used to constrain whole groups 
of theoretical scenarios. This strategy clearly generates added value to the experimental efforts, 
while at the same time it allows to achieve valuable physics results as early as possible. 


Neutrino physics led to big discoveries in the past 
decades, the greatest being the observation of neutrino 
oscillations [1], which prove that neutrino masses (albeit 
tiny) must be non-zero and that neutrino flavors mix. In 
a nutshell, this means that an electron neutrino does not 
have a fixed mass but it is rather a quantum-mechanical 
superposition of several mass eigenstates. While nowa¬ 
days most oscillation parameters are known [2] and a new 
era of precision neutrino physics has started, several fun¬ 
damental questions are still unanswered. Probably the 
most important question is whether neutrinos have a Ma- 
jorana nature, i.e., if they are identical to their antiparti¬ 
cles, which would signal a violation of lepton number and 
thus lead beyond the very successful standard model of 
particle physics. Such questions can be answered by the 
observation of neutrinoless double-/? decay (Ov/30) [3], a 
nuclear transition in which two neutrons decay simulta¬ 
neously into two protons by emitting two electrons but 
no neutrino, thus changing lepton number by two units 
and possibly signaling a Majorana neutrino mass [4]. 

The experimental search for Oi^/3/3 is a very active 
field of particle and nuclear physics. Various isotopes 
for which Ovffff is energetically allowed and many de¬ 
tection techniques are pursued. Examples are: ,6 Ge 
with high purity Ge detectors, 130 Te with Te02 bolo- 
metric detectors [5], 136 Xe with liquid Xe time projec¬ 
tion chambers [6], or Xe-loaded organic liquid scintilla¬ 
tor detectors [7]. Historically, 76 Ge-based experiments 
have been leading the field, and the resulting constraints 
on the half-life of the process are among the most strin¬ 
gent ones [8-10]. Two 76 Ge-based experiments are cur¬ 
rently active and will yield results in the near future: 
Gerda [11] and Majorana [12]. These two collabo¬ 
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rations conceive of eventually realizing a common large 
scale 76 Ge (LSGe) experiment [13], capable of probing 
the theoretically allowed region for the inverted mass or¬ 
dering (i.e., the experimentally favorable scenario corre¬ 
sponding to the upper yellow band in FIG. 1). For such 
a challenging experiment, a highly modular design and 
a staged approach implementation are needed, meaning 
that the target mass will be progressively increased. 

This paper presents realistic projections of the sensi¬ 
tivity achievable by a global analysis of the data from 
current and future experiments searching for 0 u/3/3 in 
76 Ge. Among Ouffff experiments, the ones based on ' 6 Ge 
stand out because they are designed to perform quasi 
background-free measurements. Their data can hence 
be combined without limiting assumptions on the back¬ 
ground modeling. We also point out that the sensitiv¬ 
ity of a global analysis should be considered when plan¬ 
ning the mass-increasing strategy of a project, in order 
to maximize the benefit for both theory and experiment. 
Indeed, in case no signal will be observed, large classes of 
theoretical neutrino models can be excluded already by 
intermediate stages of an experiment. 

The physics observable accessible with 0 n/3/3 experi¬ 
ments is the effective Majorana neutrino mass \m ee \ = 
\m lC f 2 cf 3 + rn 2 sf 2 cf 3 e ta21 + m 3 sf 3 e l( - a31 ~ 2S ^\, which de¬ 
pends on sines (s) and cosines (c) of the leptonic mixing 
angles 0 l:j , the mass eigenvalues, and the phases [14]. It 
is related to the Ojz/ 3/3 half-life by [15]: 

1/T°J 2 = G 0 „\Mo„\ 2 \m ee \ 2 , (1) 

where Gq u = 2.42-10~ 26 yr -1 eV~ 2 is a phase-space fac¬ 
tor and A4o is the dimensionless nuclear matrix element 
(NME) which parametrizes the nuclear physics involved. 
The allowed range for |m ee | as a function of the small¬ 
est neutrino mass m is constrained by the experimen¬ 
tal measurements of the neutrino mixing parameters, see 
FIG. 1. Nevertheless, information about the absolute 
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FIG. 1. Most general regions for the effective mass 
|m ee |, calculated using global fit parameters (v2.0) from 
nu-fit.org [2], along with four specific values of sensitivity 
considered in the text. The broadening of the bands is due to 
nuclear physics uncertainties. The disfavored regions are the 
most optimistic bounds from 136 Xe-based experiments [6, 7] 
and Planck [16], the latter converted to the smallest neutrino 
mass and averaged between both mass orderings. 


neutrino mass that is inferred by combining all exper¬ 
imental information are affected by systematic uncer¬ 
tainties of the analysis procedure [17], the NMEs [18], 
and the mixing parameters [14]. Consequently, even pin¬ 
ning down the neutrino mass ordering - whether normal, 
mi < 7712 < m 3 (blue), or inverted, m 3 < mi < m 2 
(yellow) - is challenging. 

This situation could drastically change with additional 
input from neutrino physics. The smallness of neu¬ 
trino masses can be theoretically explained by suppres¬ 
sion mechanisms at tree- [19] or loop-level [20] and the 
large mixing angles by flavor models based on discrete 
symmetries [ 21 ], which explain them by relating their 
values to properties of finite symmetry groups. While 
many of those models yield similar predictions for the 
accessible observables - so that their experimental dis¬ 
tinction is prevented unless the precision is increased by 
about two orders of magnitude - certain classes of mod¬ 
els predict clear correlations between observables. Prime 
examples are neutrino mass sum rules [22, 23], such as 
rhi + to 2 = m 3 or 1 /mi + I/TO 3 = 2/m 2 , which corre¬ 
late the complex neutrino mass eigenvalues rhi. These 
rules are complex equations and thus deliver two pieces 
of information: a constraint on the mass scale m and 
some relation between the Majorana phases 0:21,31- The 
most extensive study available [23] investigated more 
than 50 flavor models divided into 12 classes, which - 
as FIG. 2 shows - can greatly decrease the allowed range 
for |m ee |, thereby offering the possibilities of gaining valu- 
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FIG. 2. Range of the 0 u/3/3 effective mass allowed for differ¬ 
ent classes of neutrino flavor models that are characterized 
by respective sum rules. The sensitivity for different stages 
of 76 Ge-based experiments considered in the text is also dis¬ 
played. 


able knowledge on the neutrino sector already by the in¬ 
termediate steps in a staged approach towards detecting 
Oiz/3/3. 

The advantages of using high purity Ge (HPGe) detec¬ 
tors for the Oiz/3/3 search have been recognized early [24], 
HPGe detectors can be produced from germanium iso- 
topically enriched in 76 Ge ( enr Ge, typically 87% enrich¬ 
ment). The experimental signature expected for Oiz/3/3 in¬ 
side the detector is a peak in the energy spectrum at the 
Q-value of the 76 Ge decay {Qpp = 2039.061(7) keV [25]). 
Remarkable advantages of this detection technique are 
the intrinsic radio-purity of the detectors, the excellent 
spectroscopic performance (< 0 . 1 % energy resolution at 
Qpp), and the high detection efficiency. In addition, 
these detectors are a well consolidated technology widely 
used for 7 -ray spectroscopy, which proved to be reli¬ 
able and suitable for long-term experiments. The detec¬ 
tor geometries considered for (V/3/3 experiments include 
three types: coaxial, Broad Energy Germanium (BEGe), 
and P-type Point Contact (PPC) [26, 29]. Each geome¬ 
try results in a specific electric field inside the detector, 
which affects the performance of event-reconstruction 
techniques based on the time evolution of the read-out 
electrical signals (i.e., pulse shape analysis). HPGe de¬ 
tectors must be operated at cryogenic temperatures and 
are commonly installed in vacuum cryostats. This ap- 
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proach was adopted also by past Ovfd/3 experiments [8, 9], 
which operated coaxial-type detectors in ultra-low back¬ 
ground cryostats surrounded by massive lead and copper 
shieldings. 

Nowadays, the Majorana collaboration is pursuing a 
design based on PPC-type detectors and multiple cryo¬ 
stat modules built from ultrapure electroformed copper. 
Two modules are currently being assembled (i.e., the 
Majorana Demonstrator [12]) at the Sanford Under¬ 
ground Research Facility in Lead, South Dakota (USA). 
The first module hosts 16.8 kg of enr Ge detectors and will 
be fully operational in the second half of 2015. The com¬ 
pletion of the second cryostat containing further 12.6 kg 
of enr Ge detectors is scheduled by the end of 2015. The 
experiment is designed to operate the detectors at a back¬ 
ground level of 0.75-10 -3 cts/(keV- kg- yr) at Qpp. 1 The 
Gerda collaboration is instead exploring an alternative 
design in which an array of bare enr Ge detectors is oper¬ 
ated directly in ultra radio-pure liquid argon, which acts 
as coolant material, passive shielding against the exter¬ 
nal radioactivity, and active veto-system when its scin¬ 
tillation light is detected. The setup is installed in the 
Gran Sasso underground laboratories of INFN in Italy. 
Gerda has recently completed its first phase of opera¬ 
tion (Phase I), during which ~15 kg of enr Ge detectors 
(mostly of coaxial type) have been operated with a back¬ 
ground level of 10~ 2 cts/(keV- kg- yr), yielding a limit of 
T°J 2 > 2.1- 10 25 yr [10]. The apparatus is currently being 
upgraded to operate additional 17 kg of enr Ge BEGe-type 
detectors and new sensors for the argon scintillation light. 
A second data taking phase (Phase II) is planed to start 
in the second half of 2015 with a background level of 
10~ 3 cts/(keV- kg- yr) at Qpp [32]. 

The Majorana Demonstrator and Gerda 
Phase II will together start the exploration of T°7 2 at 

the scale of 10 26 yr, i.e., |m ee | ~ 0.1 eV. The results 
collected by the two experiments during the first years 
of operation will be essential to define the design of the 
LSGe experiment and down-select the best technologies 
to operate > 1000 kg of target mass at a background 
level of < 1CT 4 cts/(keV- kg- yr) at Qpp- With such 
parameters, the LSGe experiment will probe 

sensitivity at the level of 10 27 -10 28 yr and hence explore 
an essential part of the parameter space allowed for 
inverted mass ordering or with a fortunate value of 
012 and better precision on that parameter coming from 
experiments like JUNO [27] or RENO-50 [28] - even the 
whole parameter space. 

The sensitivity achievable by a global analysis of 
Gerda Phase I and Phase II, the Majorana Demon¬ 
strator, and a future LSGe experiment has been stud¬ 
ied by assuming the data sets listed in TABLE I. Follow¬ 
ing the analysis approach adopted by the Gerda collab¬ 
oration, data from Gerda Phase I are divided into two 


1 The design goal of the Majorana Demonstrator is typically 
quoted as 3 cts/(ton-yr) in a region of interest of 4keV. 


TABLE I. Parameters assumed for each data set: detector 
mass, efficiency e, background level at Qpp, energy resolution 
(full width at half maximum, FWHM) at Qpp, start time, and 
duration of the data taking At. The start time of the current 
(future) experiments is indicated with to (t 1 ) and expected to 
be in the second half of 2015 (in the 2020s). 


data 

mass 

e 

background 

FWHM 

start 

At 

set 

[kg] 


cts 

[keV] 

time 

[yr] 




keV- kg- yr 




Gerda Phase I: 






coaxial 

12.2 

0.62 

1.1- 10 -2 

4.4 

Nov 2011 

1.3 

BEGe 

2.8 

0.66 

0.5- 10" 2 

2.9 

Jul2012 

0.8 

Gerda Phase 11 






coaxial 

17.7 

0.62 

1- 10 -3 

4.0 

to 

4 

BEGe 

20.0 

0.65 

1- 10' 3 

2.5 

to 

4 

Majorana Demonstrator: 




modi 

16.8 

0.65 

0.8- 10" 3 

3.0 

to 

4 

mod2 

12.6 

0.65 

0.8-10" 3 

3.0 

fo+0.5 yr 

4 

Future large scale (LSGe) experiment: 



modi 

200 

0.65 

1- 10~ 4 

2.5 

ti 

10 

mod2 

200 

0.65 

1- 10~ 4 

2.5 

fi+lyr 

9 

mods 

200 

0.65 

1- 10~ 4 

2.5 

£i+2 yr 

8 

mod4 

200 

0.65 

1- 10~ 4 

2.5 

£i+3yr 

7 

mod5 

200 

0.65 

1- 10~ 4 

2.5 

fi+4yr 
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data sets according to the two types of detectors oper¬ 
ated. 2 The separation into two data sets is assumed also 
for Phase II. The experimental parameters such as ef¬ 
ficiencies, background level and duration are taken from 
the published values [10, 30]. The energy resolution is 
taken from the most recent R&D results [30-32], BEGe- 
type detectors will provide higher energy resolution and 
superior background reduction performance with respect 
to the coaxial type. Data from Majorana Demonstra¬ 
tor are also split between the two modules into two data 
sets. Efficiencies of PPC- and BEGe-type detectors are 
assumed to be equal. This assumption is fully consis¬ 
tent with the first results presented by the Majorana 
collaboration [26], and from which the energy resolution 
is taken. A staged approach is assumed for the LSGe 
experiments. Given realistic constraints on the produc¬ 
tion of enr Ge material 3 the total target mass of 1000 kg is 
assumed to be progressively increased by installing one 
new module with 200 kg of detectors per year. The de- 


2 The data set correspond to the “golden” and “BEGe” data sets 
of Ref. [10]. A third data set considered in the analysis of the 
collaboration (the “silver” data sets, about 6% of the overall ex¬ 
posure) is not considered here, due to its negligible contribution 
to the overall sensitivity of the experiment. 

2 The Svetlana Department facility can currently deliver 80-100 kg 

of 76 Ge per year [31]. We realistically assume that the production 

can be doubled in the next years and with some investments 
coming from the LSGe experiment. 
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tectors are considered to perform similarly to BEGe-type 
detectors. 

The total number of Oz^/3/3 events in each data set as a 
function of T °" 2 is given by: 

TV 0 " = In 2 • N a • e • r)/(m a • T 0 /,), (2) 

where Na is the Avogadro’s number, e the efficiency, r) 
the exposure, and m a the molar mass of enr Ge. In this 
work, the exposure ?y is defined as the product of total 
detector mass and data taking time. The efficiency e is 
given by the product of four contributions: the fraction 
of ,6 Ge in the detectors material (~87%), the fraction of 
the detector volume which is active (87% for coaxial, 92% 
for BEGe/PPC detectors), the efficiency of the analysis 
cuts (90%, dominated by pulse shape analysis cuts), and 
the probability that Ou/3/3 events in the detector active 
volume are correctly reconstructed at energy Qpp (92% 
for coaxial, 90% for BEGe/PPC detectors). These effi¬ 
ciencies are taken from Ref. [10]. A duty cycle of 95% 
is assumed for all experiments, which accounts for the 
time needed to calibrate the detectors and for ordinary 
hardware maintenance. 

A statistical approach is adopted to estimate the T °" 2 
lower limit achievable by a global analysis of the various 
data sets. More than 10 6 time-stamps are randomly se¬ 
lected. Given a time-stamp, background events with a 
uniform energy distribution in the range Q/3/3 A 0.1 MeV 
are generated with Monte Carlo techniques. Events are 
generated independently for each data set according to 
its background level, exposure, and efficiency. A simul¬ 
taneous fit of all data sets is hence performed, using a 
constant probability density function for the background 
and a Gaussian function for the (ffi/3/3 signal (with cen¬ 
troid at Q/ 3 / 3 ). The 90% C.L. upper limit on number of 
0 v/3(3 counts extracted from the fit is converted into a 
90% C.L. lower limit on T °" 2 by using Eq. (2). 

The fit procedure is based on an unbinned profile like¬ 
lihood analysis in which the number of Cffi/3/3 counts is 
bounded to positive values. The free parameters of the 
fit are the number of signal counts (the parameter of in¬ 
terest) and the background levels (nuisance parameters). 
Systematic uncertainties (energy scale, resolution, and 
efficiency) have been studied by adding Gaussian pull 
terms in the likelihood function and found to worsen the 
limits by < 1%. To reduce the computational time, these 
systematic uncertainties have not been included in the fi¬ 
nal simulation, as their effect is small for a limit-setting 
experiment. The coverage of the method has been tested 
for a sample of time-stamps and found to provide a con¬ 
servative overcoverage. 

The results of these computations are shown in 
FIG. 3. The top panel illustrates the integrated ex¬ 
posure over time. The increase of exposure is driven 
by Gerda Phase I (between 2012 and mid 2013), 
Gerda Phase II and the Majorana Demonstrator 
(between t 0 and to+4 yr) and the LSGe experiment (be¬ 
tween t\ and ti+10yr). The middle panel shows the 
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FIG. 3. Top panel: integrated exposure assumed for the cal¬ 
culation as a function of time before (solid line) and after 
(dashed line) efficiency correction. Middle (bottom) panel: 
distribution of the 90% C.L. lower limits on (upper limits 
on |m ee |) derived by a global analysis of multiple realizations 
of the experiments. The vertical red lines corresponds to spe¬ 
cific values of exposure/sensitivity discussed in the text. The 
time axis is broken and future dates are given with respect 
to the start of Gerda Phase If and the Majorana Demon¬ 
strator (to), and of the future LSGe experiment (fi). 

distribution of the 90% C.L. lower limits on T°" 2 . The 
uppermost part of the distribution is populated by the 
data set realizations with no background events at < 5 / 3/3 
(i.e., fully background-free). It grows linearly with the 
exposure and has a sharp cut-off due to the constraint 
Ai°" > 0 imposed on the fit. The bottom panel shows 
the distribution of the 90% upper limits on the effective 
mass \m ee \. This distribution is computed by converting 
each T°7 2 limit through Eq. (1), which introduces an ad¬ 
ditional systematic uncertainty on each limit due to the 
uncertain NME calculations. The effect of this system¬ 
atic uncertainty is maximally included in the plot assum¬ 
ing NME values in the range between 4.6 and 5.8 [23]. 4 


4 The range covers all the calculations available in the literature 
except for the shell model, which predicts an outlier value for 
NME of 2.3 [23]. Expanding the range up to the shell model 
prediction increases our upper bounds on |m ee | by a factor of 2. 
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Thus, the median line becomes a band and the central 
intervals broaden. 

FIG. 3 shows how successive experiments can be de¬ 
signed to improve the median experimental sensitivity by 
an order of magnitude. Gerda Phase I reached a sensi¬ 
tivity of Ty 2 > 2-6- 10 25 yr (|m ee | > 215-272meV) with 
an exposure of ~20 kg- yr (corresponding to the leftmost 
red line in FIG. 3). Gerda Phase II and the Majo- 
rana Demonstrator will improve the sensitivity up 
to Ty ,j > 4-10 26 yr (\m ee \ > 58-74meV) by collecting an 
exposure of 3-10 1 2 kg- yr in 4yr. The LSGe experiment 
will finally rise the sensitivity up to > 8-10 27 yr 

(|m ee | > 13-16meV) in lOyr of data taking and a final 
exposure of 8-10 3 4 5 6 7 8 9 10 kg- yr. It is noteworthy that a sensi¬ 
tivity of T°J 2 > 3-10 27 yr {\m ee \ > 19-24meV) can be 

reached with about 2-10 3 kg- yr in 4 years of data taking 
with the LSGe experiment. 

The ultimate question to answer is what can be learned 
about neutrino physics from future 76 Ge-based experi¬ 
ments. FIG. 1 shows how challenging will be to fully 
probe the parameter space allowed for \m ee \ in the most 
general situation, even considering the future LSGe and 
inverted mass ordering. Intermediate sensitivity stages 
seem not to be able to provide remarkable physics results 
unless a positive signal is observed. However, we demon¬ 
strate in FIG. 2 that whole groups of more specific neu¬ 
trino flavor models, namely those which predict particu¬ 
lar mass sum rules, can be excluded already by intermedi¬ 
ate stages. For example, the sum rule m]” 1 + rh^ 1 = mi]" 1 
yields for inverted ordering a smallest allowed neutrino 
mass of 51 meV (48 meV) for the best-fit (3<r) values 
of the neutrino mass squares. This region can be al¬ 
most probed by Gerda Phase II and the Majorana 
Demonstrator, and fully probed (i.e., even with all 
uncertainties) by first stages of the LSGe experiment. 
Thus, by using the sum rule predictions as orientation 


when planning the stages, one can exploit the synergies 
between model predictions and experimental sensitivities 
to greatly enhance the physics outcome even of the in¬ 
termediate stages. This synergy goes so far that some 
groups of models could be distinguished in spite of the 
uncertainties involved, and our considerations would be 
strengthened further by a better knowledge on the NMEs, 
on the neutrino mass ordering, or on the mixing angle 9 12 
- and even more by the observation that the sum rule 
predictions are quite stable in what regards certain types 
of theoretical (radiative) corrections [33]. Additionally, 
we would like to point out the a remarkable number of 
models could be already ruled out with ~ 2-10 3 kg- yr of 
exposure. 

Such an exposure could be collected by a single 
module of the LSGe experiment or by upgrades of 
Gerda Phase II and the Majorana Demonstrator 
which are already under consideration within the exper¬ 
imental community [13, 34]. 

In conclusion, realistic sensitivity projections have 
been presented for the current and future 76 Ge-based 
experiments. A global analysis of different data sets is 
reliable and should be performed. The global sensitiv¬ 
ity and its impact on flavor models should be carefully 
considered when designing the mass-increasing strategy 
of the future projects. 

Synergies between theory and experiment can push us 
to new frontiers in neutrino physics, provided that we 
make proper use of them. 
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